UNITUE-THEP-00/08, FAU-TP3-00/07 
http://xxx.lanl.gov/abs/nucl-th/0006082 



Nucleon Properties in the Covariant 
Quark-Diquark Model 

M. Oettel 1 , R. Alkofer 1 , and L. von Smekal 2 

1 Universitat Tubingen, Institut fur Theoretische Physik, Auf der Morgenstelle 14, 72076 Tubingen, Germany 

2 Universitat Erlangen-Nurnberg, Institut fur Theoretische Physik III, Staudtstr. 7, 91058 Erlangen, Germany 

(30 June 2000) 

Abstract. In the covariant quark-diquark model the effective Bethe-Salpeter (BS) equations for the nucleon 
and the A are solved including scalar and axialvector diquark correlations. Their quark substructure is 
effectively taken into account in both, the interaction kernel of the BS equations and the currents employed 
to calculate nucleon observables. Electromagnetic current conservation is maintained. The electric form 
factors of proton and neutron match the data. Their magnetic moments improve considerably by including 
axialvector diquarks and photon induced scalar-axialvector transitions. The isoscalar magnetic moment 
can be reproduced, the isovector contribution is about 15% too small. The ratio /j,Ge/Gm and the axial 
and strong couplings qa, g-wNN, provide an upper bound on the relative importance of axialvector diquarks 
confirming that scalar diquarks nevertheless describe the dominant 2-quark correlations inside nucleons. 
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1 Introduction 

High-precision data on nucleon properties in the medium-energy 
range are available by now or will be in near future. This is es- 
pecially the case for their electromagnetic form factors. From 
a theoretical point of view the behavior of the form factors 
indicates the necessity of a relativistic description of the nu- 
cleon. Non-relativistic constituent models generally fail beyond 
a momentum transfer of a few hundred MeV. 

From relativistic quantum mechanics of three constituent 
quarks models employing effective Hamiltonian descriptions 
were deduced. However, the necessity for the effective Hamil- 
tonian to comply with the Poincare algebra and, at the same 
time, with the covariance properties of the wave functions leads 
to fairly complicated constraints to ensure the covariance of 
current operators [1]. While some phenomenological studies 
relax these constraints and allow for covariance violations of 
one-body currents [2], others put emphasis on the consistent 
transformation of those components of one-body currents that 
are relevant in light- front Hamiltonian dynamics [3,4,5]. The 
inclusion of two-body currents was considered within a semi- 
relativistic chiral quark model in the study of ref. [6]. Generally, 
in order to describe the phenomenological dipole shape of the 
electric form factor of the proton, additional form factors for 
the constituent quarks need to be introduced in a phenomeno- 
logical way in these quantum mechanical models [5,6]. 

Models based on the quantum field theoretic bound state 
equations for baryons on the other hand have proved capable of 
describing the electric form factors of both nucleons quite suc- 
cessfully without such additional assumptions [7,8,9]. Param- 
eterizations of covariant Faddeev amplitudes of the nucleons 
for calculating various form factors were explored in refs. [8] 
employing impulse approximate currents, the field theoretic 
analog of using one-body currents. While the covariance of the 
corresponding nucleon amplitudes is of course manifest in the 
quantum field theoretic models, current conservation requires 
one to go beyond the impulse approximation, however, also 
in these studies [9,10,11]. Furthermore, the invariance under 
(4-dimensional) translations ramifies into certain properties of 
the baryonic bound state amplitudes which are generally not 
reflected by the parameterizations but result only for solutions 
to their quantum field theoretic bound state equations such as 
those obtained in ref. [9] . 

The axial structure of the nucleon is known far less pre- 
cisely than the electromagnetic one. The theoretical studies 
have mainly been focused on the soft point limit. Even though 
precise experimental data on the pion-nucleon and the axial 
form factor for finite Q 2 are difficult to obtain, and thus prac- 
tically unavailable, it would be very helpful to compare the 
various theoretical results (see, e.g., refs. [7] and [8]) to such 
data. 

In this paper, we investigate the structure of the nucleon 
within the covariant quark-diquark model. In previous applica- 
tions of this model, including calculations of quark distribution 
functions [12] and of various nucleon form factors in the im- 
pulse approximation [7], pointlike scalar diquark correlations 
were employed. The octet and decuplet baryon spectrum is 
described in ref. [13] maintaining pointlike scalar and axialvec- 
tor diquarks. A correct description of the nucleons' electric 
form factors and radii was obtained in ref. [9] by introduc- 
ing diquark substructure. We extend this latter calculation to 
non-pointlike axialvector diquarks which allows us to include 



the Zi-resonance in the calculations. The model is fully rela- 
tivistic, reflects gauge invariance in the presence of an elec- 
tromagnetic field, and it allows a direct comparison with the 
semi-relativistic treatment obtained from the Salpeter approx- 
imation (to the relativistic Bethe-Salpeter equation employed 
herein). This comparison was done in ref. [14] showing that 
both, the static observables (magnetic moments, pion-nucleon 
coupling constant and axial coupling constant) and the be- 
haviour of the neutron electric form factor (or its charge ra- 
dius, which is particularly sensitive to the specific assumptions 
of any nucleon model), in such a treatment deviate consid- 
erably from the fully relativistic one. These results indicate 
that further studies of the quantum field theoretic bound state 
equations for nucleons are worth pursuing. 

In the next section, we briefly review the basic notions of 
the quark-diquark model. In Section 3 the model expressions 
for the electromagnetic, the pion-nucleon and the weak-axial 
form factors are derived. Starting from the quark vertices and 
their respective currents we construct effective diquark vertices 
and fix their strengths by resolving the quark-loop structure 
of the diquarks. We employ two parameter sets in our calcula- 
tions. One set is chosen to fit nucleon properties alone whereas 
the other one includes the mass of the delta resonance. While 
the electric form factors are essentially identical for both sets, 
differences occur in the magnetic form factors. We present our 
results and discuss their implications in section 4. The electro- 
magnetic form factors are thereby compared to experimental 
data. In conclusion we comment on future perspectives within 
this framework. 



2 The Quark-Diquark Model 

Nucleon and delta are modelled as bound states of three con- 
stituent quarks. In order to make the relativistic three-body 
problem tractable, we neglect any 3-particle irreducible inter- 
actions between the quarks and assume separable correlations 
in the two-quark channel. The latter assumption introduces 
non-pointlike diquark correlations. The first assumption allows 
to derive a relativistic Faddeev equation for the 6-point quark 
function and the assumed separability reduces it to an effective 
quark-diquark Bethe-Salpeter (BS) equation. In the following, 
we work in Euclidean space. 

In this article we restrict ourselves to scalar and axialvector 
diquarks which are introduced, as stated above, via the sep- 
arability assumption for the connected and truncated 4-point 
quark function: 

xUp)d^(p)xMi)- (2-1) 

P is the total momentum of the incoming and the outgoing 
quark pair, p and q are the relative momenta between the 
quarks in the two channels. The propagators of scalar and axi- 
alvector diquark in eq. (2.1) are those of free spin-0 and spin-1 
particles, 

*w - -whr, ■ < 2 - 2 ' 

Correspondingly, Xa&{p) an d Xapip) are the respective quark- 
diquark vertex functions. Here, we maintain only their dom- 
inant Dirac structures which are multiplied by an invariant 
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function P(p 2 ) of the relative momentum p between the quarks 
to parameterize the quark substructure of the diquarks. The 
Pauli principle then fixes this relative momentum to be an- 
tisymmetric in the quark momenta p a and pp [9], i.e., p = 
\{Pa — Pp). Besides the structure of the vertex functions in 
Dirac space, they belong to the anti-triplet representation in 
color space, i.e. they are proportional to cabd with color in- 
dices A, B for the quarks and D labelling the color of the di- 
quark. Furthermore, the scalar diquark is an antisymmetric fla- 
vor singlet represented by (r2) a &, and the axialvector diquark is 
a symmetric flavor triplet which can be represented by (T2Tk) a b- 
Here, a and 6 label the quark flavors and k the flavor of the 
axialvector diquark. Thus, with all these indices made explicit, 
the vertex functions read 1 



I \ I D / •> ( T '2)ab £ABD 

X<x0(p) = M7 C%/3 P(p) 



XUP) = 9a(^CU P(P) 



V2 V2 ' 

(T~2Tk)ab eABD 
y/2 ^/2 



(2.4) 
(2.5) 



Here, C denotes the charge conjugation matrix; g a and g s are 
effective coupling constants of two quarks to scalar and axi- 
alvector diquarks, respectively. 

For the scalar function P(p), we employ a simple dipole 
form with an effective width A which has proven very successful 
in describing the phenomenological dipole form of the electric 
form factor of the proton in ref. [9] , 



P(p) 



A 



A 2 + p' 2 



(2.6) 



This models the non-pointlike nature of the diquarks. It fur- 
thermore provides for the natural ultraviolet regularity of the 
interaction kernel in the nucleon and delta BS equations to be 
derived below. 

We can compute the coupling constants g s and g a by put- 
ting the diquarks on-shell and evaluating the canonical nor- 
malization condition by using that the vertex functions Xafiip) 
and Xapip) can De viewed as diquark amplitudes with trun- 
cated quark legs: 



I f d *P - 
4 J (2tt) 4 Xa(: 



d 



(o-) 



Xis = 2m s 



d 4 p 
(2tt)« 



(Xa/3 



9PM <*7,/3; 



(2.7) 
(2.8) 



Hereby G^^ s is the free propagator for two quarks, sym- 
metrized (+) for the axialvector diquark and antisymmetrized 
(— ) for the scalar diquark [9]. Furthermore, {x v ) T = x" ~ 
P v (P fl X fl ) is the transverse part of the vertex function x" 
(note that the pole contribution of the axialvector diquark is 
transverse to its total momentum, and the sum over the three 
polarization states provides an extra factor of 3 on the r.h.s. 
of eq. (2.8). as compared to eq. (2.7)). With normalizations 
as chosen in eqs. (2.4,2.5) the traces over the color and flavor 
parts yield no additional factors. 



1 Symbolically denoting the totality of quark indices by the 
same Greek letters that are used as their Dirac indices should 
not create confusion. The particular flavor structures are tied 
to the Dirac decomposition of the diquarks, color-3 is fixed. 



2.1 Nucleon Amplitudes 

The nucleon BS amplitudes (or wave functions) can be de- 
scribed by an effective multi-spinor characterizing the scalar 
and axialvector correlations, 



V(p,P)u(P,a) 



V 5 (p,P) 



u(P,s). 



(2.9) 



u(P,s) is a positive-energy Dirac spinor (of spin s), p and P 
are the relative and total momenta of the quark-diquark pair, 
respectively. The vertex functions are defined by truncation of 
the legs, 



D- 1 








(2.10) 



The diquark propagators D and D^ v are defined in eqs. (2.2,2.3) 
and S _1 denotes the inverse quark propagator in eq. (2.10). In 
the present study we employ that of a free constituent quark 
with mass m„, 



S~\p) = -if 



(2.11) 



The coupled system of BS equations for the nucleon ampli- 
tudes or their vertex functions can be written in the following 
compact form, 



f d 4 p' x , (9* 



(p',P)=0, 



(2.12) 



in which G^ 1 (p,p' , P) is the inverse of the full quark-diquark 
4-point function. It is the sum of the disconnected part and 
the interaction kernel. 

Here, the interaction kernel results from the reduction of 
the Faddeev equation for separable 2-quark correlations. It de- 
scribes the exchange of the quark with one of those in the 
diquark which is necessary to implement Pauli's principle in 
the baryon. Thus, 

G-\p,p',P) = 



(2tt) 4 5 4 {p-p')S-\p q ) 



D-\p d ) 



(£>" »r\Pd) 

-X(p1) S T (q) X(pT) V3 X M '(p1) S T (q) X(p1) 

V3 X (pD S T ( q ) x m (p?) X m '(pI) S T ( q) x"(p?) 



(2.13) 



Herein, the flavor and color factors have been taken into ac- 
count explicitly, and x, X M stand for the Dirac structures of 
the diquark-quark vertices (multiplied by the invariant func- 
tion P(pi r 2 j), see ec l- (2-6)). The freedom to partition the total 
momentum between quark and diquark introduces the param- 
eter r) G [0, 1] with p q = rjP + p and p d = (1 — n)P — p. 
The momentum of the exchanged quark is then given by q — 
—p — p' + (1 — 27/) P. The relative momenta of the quarks in 
the diquark vertices x and x are p2 = p + p' /2 — (1 — 3ri)P/2 
and pi — p/2 + p' — (1 — 3r/)P/2, respectively. Invariance un- 
der (4- dimensional) translations implies that for every solution 
$ (p, P; 771 ) of the BS equation there exists a family of solutions 
of the form <P(p + (772 — iJi)P, P', ^2). 

Using the positive energy projector with nucleon bound- 
state mass M„, 



2 



1 + 



iM n 



(2.14) 
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the vertex functions can be decomposed into their most general 
Dirac structures, 



(2.15) 



(2.16) 



iM n 



In the rest frame of the nucleon, P = (0,iM„), the unknown 
scalar functions Si and Ai are functions of p 2 = p M p M and of 
the angle variable z — P ■ p, the cosine of the (4-dimensional) 
azimuthal angle of p M . Certain linear combinations of these 
eight covariant components then lead to a full partial wave 
decomposition, see ref. [13] for more details and for examples 
of decomposed amplitudes assuming pointlike diquarks. These 
nucleon amplitudes have in general a considerably broader ex- 
tension in momentum space than those obtained herein with 
including the quark substructure of diquarks, however. 

The BS solutions are normalized by the canonical condition 

d 



(2k) 
*(p',P, 



pfi 9^G-\ P ',p,P) 



*(p,Pn) 



2.2 Delta Amplitudes 

The effective multi-spinor for the delta baryon representing 
the BS wave function can be characterized as ^ A (p,P)u" '(P) 
where u v (P) is a Rarita-Schwinger spinor. The momenta are 
defined analogously to the nucleon case. As the delta state is 
flavor symmetric, only the axialvector diquark contributes and, 
accordingly, the corresponding BS equation reads, 



j4 / 
a p 



(2n) 



1 G A 1 (p,p',P)W^(p',P)^0 



(2.18) 
in the A- 



where the inverse quark-diquark propagator G A 
channel is given by 

G A \p,p',P) = (2^) 4 5 4 ( P -p') S-'M (D^')-\ Pd ) + 

X»\pl)S T {q)r{pX)- (2.19) 

The general decomposition of the corresponding vertex func- 
tion <P A V , obtained as in eq. (2.10) by truncating the quark and 
diquark legs of the BS wave function W^f , reads 



iD, + w -JD 2 )A^ 



iM A 



(£i 



M A iM A 



iM A y M A F b 'iM A 



A 



(2.20) 



Here, A^ v is the Rarita-Schwinger projector, 

i_P^j_ 
3 



A" v = A+ 



1 M v 2P»P» 
3 7 7 3 M 2 



P"Y 



M A 



which obeys the constraints 
P"/!"" = ^A^ = 0. 



(2.21) 



(2.22) 





Set I 


Set II 


m q 


[GeV] 


0.360 


0.425 


m sc 


[GeV] 


0.625 


0.598 


m a x 


[GeV] 


0.684 


0.831 


A 


[GeV] 


0.95 


0.53 


9s 




9.29 


22.10 


9a 




6.97 


6.37 


M A 


[GeV] 


1.007 


1.232 


M„ 


[GeV] 


0.939 


0.939 



Table 2.1. The two parameter sets employed in the calcula- 
tions herein together with the values of couplings and bound 
state masses obtained with these sets. 



Therefore, the only non-zero components arise from the con- 
traction with the transverse relative momentum p pT = p M — 
P fl (p-P). The invariant functions Di and Ei in eq. (2.20) again 
depend on p 2 and p ■ P. The partial wave decomposition in the 
rest frame is given in ref. [13], and again, the Z\-amplitudes 
from pointlike diquarks of [13] are wider in p 2 than those ob- 
tained herein. 



2.3 Solutions for BS Amplitudes of Nucleons and A 

The nucleon and A BS equations are solved in the baryon 
rest frame by expanding the unknown scalar functions in terms 
of Chebyshev polynomials [13,9]. Iterating the integral equa- 
tions yields a certain eigenvalue which by readjusting the pa- 
rameters of the model is tuned to one. Altogether there are 
four parameters, the quark mass m q , the diquark masses m 3C 
and max and the diquark width A. 

In the calculations presented herein we shall illustrate the 
consequences of our present model assumptions with two dif- 
ferent parameter sets as examples which emphasize slightly dif- 
ferent aspects. For Set I, we employ a constituent quark mass 
of m q — 0.36 GeV which is close to the values commonly used 
by non- or semi-relativistic constituent quark models. Due to 
the free-particle poles in the bare quark and diquark propaga- 
tors used presently in the model, the axialvector diquark mass 
is below 0.72 GeV and the delta mass below 1.08 GeV. On 
the other hand, nucleon and delta masses are fitted by Set II, 
i.e., the parameter space is constrained by these two masses. 
In particular, this implies m q > 0.41 GeV. Both parameter 
sets together with the corresponding values resulting for the 
effective diquark couplings and baryon masses are given in ta- 
ble 2.1. 

Two differences between the two sets are important in the 
following: The strength of the axialvector correlations within 
the nucleon is rather weak for Set II, since the scalar diquark 
contributes 92% to the norm integral of eq. (2.17) while the ax- 
ialvector correlations and scalar-axialvector transition terms 
together give rise only to the remaining 8% for this set. For 
Set I, the fraction of the scalar correlations is reduced to 66%, 
the axialvector correlations are therefore expected to influence 
nucleon properties more strongly for Set I than for Set II. Sec- 
ondly, the different constituent quark masses affect the mag- 
netic moments. We recall that in non-relativistic quark mod- 
els the magnetic moment is roughly proportional to M n /m q 
and that most of these models thus employ constituent masses 
around 0.33 GeV. 
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Fig. 3.1. Impulse approximate contributions to the electro- 
magnetic current. For the scalar-axialvector transition, a dia- 
gram analogous to the third one (with initial and final nucleon 
states interchanged) has to be computed. 

3 Observables 

3.1 Electromagnetic Form Factors 

The Sachs form factors Ge and Gm can be extracted from the 
solutions of the BS equations using the relations 



Ge 



2P 2 



Tr{.r)P», Gm = ^ Tr<J">7£ 



(3.1) 



where P = (Pi + P f )/2, ^ = 7 M - P^tf , and the spin-summed 
matrix element (J M ) is given by 



(,F) = (P f ,s f \.r\P, Sl ) Y, u{P s ,s s )u{Pi, Si ) 



f d 4 Pf [ 
J (2^) 4 J 



d 4 Pi ,f, 



(2k) 



¥(Pf,Pf)J" #{Pi,Pi)- 



(3.2) 



The current J M herein is obtained as in ref . [9] . It represents a 
sum of all possible couplings of the photon to the inverse quark- 
diquark propagator G _1 given in eq. (2.13). This construction 
which ensures current conservation can be systematically de- 
rived from the general "gauging technique" employed in refs. 
[10,11]. 

The two contributions to the current that arise from cou- 
pling the photon to the disconnected part of G _1 , the first term 
in eq. (2.13), yield the impulse approximate couplings to quark 
and diquark. They are graphically represented by the middle 
and the upper diagram in figure 3.1. The corresponding ker- 
nels, to be multiplied by the charge of the respective quark or 
diquark upon insertion into the r.h.s. of eq. (3.2), read, 

r q = (27T) 4 8 4 (p f ^ Pl - V Q)r^D- 1 (k d ), (3.3) 

= ( 27 4 s '(Pf - Pi + (i - i)Q)rZ$H s-\k q ). (3.4) 

Here, the inverse diquark propagator D^ 1 comprises both, 
scalar and axialvector components. The vertices in eqs. (3.3) 



and (3.4) are the ones for a free quark, a spin-0 and a spin-1 
particle, respectively, 



1 1 



-i 7 M , /T c = -(pd + fcd)", and (3.5) 
= -(Pd + kdY 5 a0 + P a d 5^ + k$ 8» a + 

n {Q 8^ a - Q a 8^). (3.6) 

The Dirac indices a, j3 in (3.6) refer to the vector indices of the 
final and initial state wave function, respectively. The axialvec- 
tor diquark can have an anomalous magnetic moment k. We 
obtain its value from a calculation for vanishing momentum 
transfer (Q 2 = 0) in which the quark substructure of the di- 
quarks is resolved, i.e., in which a (soft) photon couples to the 
quarks within the diquarks. The corresponding contributions 
are represented by the upper and the right diagram in fig- 
ure 3.2. The calculation of k is provided in Appendix A.l. The 
values obtained from the two parameter sets are both very close 
to k = 1 (see table A.l). This might seem understandable from 
nonrelativistic intuition: the magnetic moments of two quarks 
with charges qi and q-z add up to (qi + q2)/m q , the magnetic 
moment of the axialvector diquark is (1 + n)(qi +q2)/m a x and 
if the axialvector diquark is weakly bound, m ax ^ 2m q , then 
k ±2 1. In the following we use k = 1. 

The vertices in eqs. (3.5) and (3.6) satisfy their respective 
Ward-Takahashi identities, i. e. those for free quark and diquark 
propagators (c./., eqs. (2.2,2.3) and (2.11)), and thus describe 
the minimal coupling of the photon to quark and diquark. 

We furthermore take into account impulse approximate 
contributions describing the photon-induced transitions 
between scalar and axialvector diquarks as represented by the 
lower diagram in figure 3.1. These yield purely transverse cur- 
rents and do thus not affect current conservation. The tensor 
structure of these contributions resembles that of the triangle 
anomaly. In particular, the structure of the vertex describing 
the transition from axialvector (with index 0) to scalar diquark 
is given by 



1 Sf). ° 



(3.7) 



and that for the reverse transition from an scalar to axialvector 
(index a) by, 



r'' 

- 1 n 



^sa fj,apX 



2M n 



( P d + k d ) p Q\ 



(3.8) 



The tensor structure of these anomalous diagrams (for Q — > 0) 
is derived by resolving the diquarks in Appendix A.l in a way 
as represented by the lower diagram in figure 3.2. The explicit 
factor 1/M n was introduced to isolate a dimensionless constant 
K sa . Its value is obtained roughly ~ 2.1 (with the next 

digit depending on the parameter set, c.f. table A.l). 

Upon performing the flavor algebra for the current matrix 
elements of the impulse approximation, one obtains the follow- 
ing explicit forms for proton and neutron, 



/ rftiimp " I j[J,\sc — sc . ^ i jfi \ sc — sc . / jfi \ ax— ax 
W )p - 3W9/ + gWsc) +{Jax) 

^ «jr a > sc - ax + ( j^ s ) ax - sc ) , 



(3.9) 



<.#„>"-") - ^ «jr a > sc - ax + (j^r- sc ) .(3.10) 
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~ (i,ap 
*■ ax 



*■ sa 



Fig. 3.2. Resolved vertices: photon-scalar diquark, photon- 
axialvector diquark and anomalous scalar-axialvector diquark 
transition. 




The superscript 'sc-sc' indicates that the current operator is to 
be sandwiched between scalar nucleon amplitudes for both the 
final and the initial state in eq. (3.2). Likewise 'sc-ax' denotes 
current operators that are sandwiched between scalar ampli- 
tudes in the final and axialvector amplitudes in the initial state, 
etc.. We note that the axialvector amplitudes contribute to the 
proton current only in combination with diquark current cou- 
plings. 

Current conservation requires that the photon also be cou- 
pled to the interaction kernel in the BS equation, i.e., to the 
second term in the inverse quark-diquark propagator G~ of 
eq. (2.13). The corresponding contributions were derived in [9]. 
They are represented by the diagrams in figure 3.3. In particu- 
lar, in [9] it was shown that, in addition to the photon coupling 
with the exchange-quark (with vertex r q ), irreducible (seagull) 
interactions of the photon with the diquark substructure have 
to be taken into account. The structure and functional form 
of these diquark-quark-photon vertices is constrained by Ward 
identities. The explicit construction of ref. [9] yields seagull 
couplings of the following form (with M M denoting that for the 
scalar diquark and M^^ that for the axialvector with Lorentz 
index (3): 



^f^[^( Pl+ ^)-X^( Pl )].(3.H) 

Here, q q denotes the charge of the quark with momentum p q , 
q ex the charge of the exchanged quark with momentum q 1 , and 
pi is the relative momentum of the two, pi = (p q — q')/2 (see 
figure 3.3). The conjugate vertices M 11 ''"' are obtained from 



the conjugation of the diquark amplitudes \ 



together with the replacement pi 
ref. [9]). 



P2 



in eq. (3.11) 
- k q )/2 (c./., 



Fig. 3.3. Exchange quark and seagull diagrams. 



Regarding the numerical evaluation of the diagrams we re- 
mark that in the case of the impulse approximation diagrams 
we use the covariant decomposition of the vertex function $ 
given in eqs. (2.15,2.16) together with the numerical solution 
for the scalar functions Si, Ai, Bi, d. The continuation of 
these functions from the nucleon rest frame to the Breit frame 
is described in detail in ref. [9]. For finite momentum transfer, 
care is needed in treating the singularities of the quark and di- 
quark propagators that appear in the single terms of eq. (3.2). 
In ref. [9] it was shown that for some kinematical situations 
explicit residues have to be taken into account in the calcula- 
tion of the impulse approximation diagrams. This applies also 
to the calculations presented here. 

The computation of the diagrams given in figure 3.3 in- 
volves two four-dimensional integrations. As the singularity 
structure of these diagrams becomes quite intricate, we re- 
sorted to an expansion of the wave function 9 analogously 
to the expansion of The corresponding scalar functions that 
have been computed in the rest frame as well show a much 
weaker convergence in terms of the Chebyshev polynomials 
than the scalar functions related to the expansion of the vertex 
function [9]. As a result, the numerical uncertainty (for these 
diagrams only) exceeded the level of a few percent beyond mo- 
mentum transfers of 2.5 GeV 2 . Due to this limitation (which 
could be avoided by increasing simply the used computer time) 
the form factor results presented in section 4 are restricted to 
the region of momentum transfers below 2.5 GeV 2 . 

These remarks concerning the numerical method apply also 
to the computation of the strong form factor g n NN and the 
weak form factor gA which will be discussed in the following 
subsection. 
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3.2 The Strong Form Factor q-knn and the Weak Form 
Factor qa 

The coupling of the pion to the nucleon, described by a 
pseudoscalar operator, and the pseudovector currents of weak 
processes such as the neutron /3-decay are connected to each 
other in the soft limit by the Goldberger-Treiman relation. 

The (spin-summed) matrix element of the pseudoscalar 
density Jf can be parameterized as 



(J? ) = A + (P f )r a To g^ NN (Q 2 )A+(PA. 



(3.12) 



Straightforward Dirac algebra allows to extract form factor as 
the following trace, 



<?7rJVAr(Q 2 ) 



2M 2 r a 
-^Tr<J 5 > 75T . 



(3.13) 



To compute the form factors we first specify a suitable quark- 
pion vertex, evaluate an impulse approximate contribution cor- 
responding to the upper diagram of figure 3.1, and an exchange 
contribution analogous to the upper diagram of figure 3.3. The 
structures and strengths (for Q — > 0) for the couplings of the 
diquarks to the pion and the axialvector current (the remain- 
ing two impulse approximation diagrams in figure 3.1) are ob- 
tained from resolving the diquarks in a way similar to their 
electromagnetic couplings in Appendix A. 2. 

The structure of the inverse quark propagator, given by 
S~ 1 {p) = — ' l i> A(p 2 ) — B(p 2 ), suggests that we use for the pion- 
quark vertex 



r 5 = -75 -7-t 



(3.14) 



and discard the three additionally possible Dirac structures (f n 
is the pion decay constant). The reason is that in the chiral 
limit eq. (3.14) is the exact pion BS amplitude for equal quark 
and antiquark momenta, since the Dyson- Schwinger equation 
for the scalar function B agrees with the BS equation for a 
pion of zero momentum in this limit. Of course, the subdom- 
inant amplitudes should in principle be included for physical 
pions (with momentum —P 2 — m\), when solving the Dyson- 
Schwinger equation for A and B and the BS equation for the 
pion in mutually consistent truncations [15,16,17]. Herein we 
employ A(p 2 ) — 1 and B(p 2 ) = m q . 

The matrix elements of the pseudovector current are para- 
meterized by the form factor qa{Q 2 ) and the induced pseu- 
doscalar form factor gp{Q 2 ), 

= A + {P t ) £ [i^To9A{Q 2 ) + Q^ 5 g P (Q 2 )] A + (P). 

(3.15) 

For Q 2 — > the Goldberger-Treiman relation, 

9a(0) = f*g*NN(0)/M n , (3.16) 

then follows from current conservation and the observation 
that only the induced pseudoscalar form factor gp(Q 2 ) has 
a pole on the pion mass-shell. 

By definition, gA describes the regular part of the pseu- 
dovector current and gp the induced pseudoscalar form factor. 
They can be extracted from eq. (3.15) as follows: 

2M n 



9a(Q 2 ) = 



757 - 175- 



-Q' 



gp(Q 2 ) = 



2M n 

Q 2 



9a(Q 2 ) 



Mr, 

Q 2 



Tr(J 5 a ^)Q' i 75 r a 



(3.17) 
(3.18) 



We again use chiral symmetry constraints to construct the 
pseudovector-quark vertex. In the chiral limit, the Ward-Ta- 
kahashi identity for this vertex reads, 



£(S- 1 (fc) 7 5 +75S- 1 (p)) , (Q = k-p). 



(3.19) 



To satisfy this constraint we use the form of the vertex pro- 
posed in ref. [18], 



rr = -n M 75- 



+ 



(3.20) 



The second term which contains the massless pion pole does 
not contribute to gA as can be seen from eq. (3.17). From these 
quark contributions to the pion coupling and the pseudovector 
current alone, eqs. (3.20) and (3.13) would thus yield, 



( ,; ^(Q 2 ) = ^^(0). 



lim 



(3.21) 



Here, the Goldberger-Treiman relation followed if the pseu- 
dovector current was conserved or, off the chiral limit, from 
PCAC. 

Current conservation is a non-trivial requirement in the rel- 
ativistic bound state description of nucleons, however. First, we 
ignored the pion and pseudovector couplings to the diquarks 
in the simple argument above. For scalar diquarks alone which 
themselves do not couple to either of the two, 2 pseudovector 
current conservation could in principle be maintained by in- 
cluding the couplings to the interaction kernel of the nucleon 
BS equation in much the same way as was done for the elec- 
tromagnetic current. 

Unfortunately, when axialvector diquarks are included, even 
this will not suffice to maintain current conservation. As ob- 
served recently in ref. [11], a doublet of axialvector and vector 
diquarks has to be introduced, in order to comply with chi- 
ral Ward identities in general. The reason essentially is that 
vector and axialvector diquarks mix under a chiral transfor- 
mation whereas this is not the case for scalar and pseudoscalar 
diquarks. Since vector diquarks on the one hand introduce six 
additional components to the nucleon wave function, but are on 
the other hand not expected to influence the binding strongly, 
here we prefer to neglect vector diquark correlations and to 
investigate the axial form factor without them. 

The pion and the pseudovector current can couple to the 
diquarks by an intermediate quark loop. As for the anoma- 
lous contributions to the electromagnetic current, we derive the 
Lorentz structure of the diquark vertices and calculate their ef- 
fective strengths from this quark substructure of the diquarks 
in Appendix A. 2. 

As mentioned above, no such couplings arise for the scalar 
diquark. The axialvector diquark and the pion couple by an 
anomalous vertex. Its Lorentz structure is similar to that for 
the photon-induced scalar-to-axialvector transition in eq. (3.7): 



U p ^ v ( Pi + k d fQ v (1 - 25 a2 )ie abc . (3.22) 



■pp\,abc ^ax ^ 

5 > a * ~ 2M~J n 

Here, a is the flavor index of the pion or, below, of the pseu- 
dovector current, while 6 and c are those of the outgoing and 
the incoming axialvector diquarks (with Lorentz indices p and 
A) according to eq. (2.5), respectively. The factor m q / f-n comes 

2 This can be inferred from parity and covariance. 
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from quark-pion vertex (3.14) in the quark- loop (see Appendix 
A. 2), and the nucleon mass was introduced to isolate a dimen- 
sionless constant K ax . 

The pseudovector current and the axialvector diquark are 
also coupled by anomalous terms. As before, we denote with p 
and A the Lorentz indices of outgoing and incoming diquark, 
respectively, and with p the pseudovector index. Out of three 
possible Lorentz structures for the regular part of the ver- 
tex, p^ {"^PaQ 13 , e" pXa Q a and e MpA "(p d + k d ) a , in the limit 
Q — > only the last term contributes to qa- We furthermore 
verified numerically that the first two terms yield negligible 
contributions to the form factor also for finite Q. Again, the 
pion pole contributes proportional to Q M , and our ansatz for 
the vertex thus reads 

rZl' abC = %^ pA >d + k d y 1(1 - 2S a2 )ie abc + 

^UnX hC - (3-23) 

For both strengths we roughly obtain K ax ~ K^ ax ~ 4.5 
slightly dependent on the parameter set, see table A. 2 (in Ap- 
pendix A. 2). 

Scalar-to-axialvector transitions are also possible by the 
pion and the pseudovector current. An effective vertex for the 
pion-mediated transition has one free Lorentz index to be con- 
tracted with the axialvector diquark. Therefore, two types of 
structures exist, one with the pion momentum Q, and the other 
with any combination of the diquark momenta p d and k d . If we 
considered this transition as being described by an interaction 
Lagrangian of scalar, axialvector and pseudoscalar fields, terms 
of the latter structure were proportional to the divergence of 
the axialvector field which is a constraint that can be set to 
zero. We therefore adopt the following form for the transition 
vertex, 

rs£ = -i£ a ^Q p (2<r 2 - i)s ab . (3.24) 

The flavor and Dirac indices of the axialvector diquark are b 
and p. This vertex corresponds to a derivative coupling of the 
pion to scalar and axialvector diquark. 

The pseudovector-induced transition vertex has two Lorentz 
indices, denoted by p for the pseudovector current and p for 
the axialvector diquark. From the momentum transfer Q p and 
one of the diquark momenta, altogether five independent ten- 
sors can be constructed: 8"", Q P Q P , Q p p p d , p^Q" and p p p p d 
(the totally antisymmetric tensor has the wrong parity). We 
assume as before, that all terms proportional to Q M are con- 
tained in the pion part (and do not contribute to qa)- From 
the diquark loop calculation in Appendix A. 2 we find that the 
terms proportional to p d Q p and p d p d can again be neglected 
with an error on the level of one per cent. Therefore, we use a 
vertex of the form, 

mf = iM nK i, m s pp i(2r 2 - i)s ab + ^unf a . 

(3.25) 

For the strengths of these two transition vertices we obtain 
Ks a ~ 3.9 and (on average) «^ j3a ~ 2.1, see table A. 2 in Ap- 
pendix A. 2. 

We conclude this section with the observation that eq. (3.21) 
remains valid with all contributions from quarks and diquarks 



included in the pseudoscalar density (3.12) and the pseudovec- 
tor current (3.15) of the nucleon. The extension of this state- 
ment to include the diquark couplings follows from the fact 
that their vertices for the pseudovector current, eqs. (3.23) 
and (3.25), contain the pion pole in the form entailed by their 
Ward identities. Eq. (3.21) is then obtained straightforwardly 
by inserting g n NN from eq. (3.13) term by term into the cor- 
responding equation (3.18) for gp and expanding to leading 
order in Q 2 . 



4 Results 

4-1 Electromagnetic Form Factors 

The results for the electric form factors are shown in figure 
4.1. The phenomenological dipole behavior of the proton Ge 
is well reproduced by both parameter sets. The electric radius 
(see table 4.1) is predominantly sensitive to the width of the BS 
wave function. This is different for the neutron. Here, stronger 
axialvector correlations tend to suppress the electric form fac- 
tor as compared to the calculation of ref. [9] where only scalar 
diquarks were maintained. A smaller binding energy compen- 
sates this effect, so that the results for Ge are basically the 
same for both parameter sets, see figure 4.1. 

For the magnetic moments summarized in table 4.2 two 
parameters are important. First, in contrast to non-relativistic 
constituent models, the dependence of the proton magnetic 
moment on the ratio M n /m q is stronger than linear. As a re- 
sult, the quark impulse contribution to p p with the scalar di- 
quark being spectator, which is the dominant one, yields about 
the same for both sets, even though the corresponding nucleon 
amplitudes of Set I contribute about 25% less to the norm 
than those of Set II. Secondly, the scalar-axialvector transi- 
tions contribute equally strong (Set I) or stronger (Set II) than 
the spin flip of the axialvector diquark itself. While for Set II 
(with weaker axialvector diquark correlations) the magnetic 
moments are about 30% too small, the stronger diquark cor- 
relations of Set I yield an isovector contribution which is only 
15% below and an isoscalar magnetic moment slightly above 
the phenomenological value. 

Stronger axialvector diquark correlations are favorable for 
larger values of the magnetic moments as expected. If the 
isoscalar magnetic moment is taken as an indication that those 
of Set I are somewhat too strong, however, a certain mismatch 
with the isovector contribution remains, also with axialvector 
diquarks included. 

Recent data from ref. [24] for the ratio p p Ge /Gm is com- 
pared to our results in figure 4.2. The ratio obtained from Set 
II with weak axialvector correlations lies above the experimen- 
tal data, and that for Set I below. The experimental observa- 





Set I 


Set II 


expt. 


(r p 


cl 


[fm] 


0.88 


0.81 


0.836 ±0.013 


(rl 


)el 


[fm 2 ] 


-0.12 


-0.10 


-0.113 ±0.007 


(r p 


mag 


[fm] 


0.84 


0.83 


0.843 ± 0.013 


(r„ 


) mag 


[fm] 


0.84 


0.83 


0.840 ± 0.042 



Table 4.1. Nucleon electric and magnetic radii for the two 
parameter sets compared to the experimental values from [19] 
(for (r 2 ) c i) and [20] (for the remaining radii). 
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Q 2 [ GeV 2 ] Q 2 [ GeV 2 ] 

Fig. 4.1. Electric form factors of both nucleons for the parameter sets of table 2.1. Experimental data for the proton is 
taken from [20]. The older neutron data analysis [21] contains more systematic uncertainties (due to specific nucleon-nucleon 
potentials) than the more recent data from [22,23]. 



tion that this ratio decreases significantly with increasing Q 2 
(about 40% from Q 2 = to 3.5 GeV 2 ), can be well reproduced 
with axialvector diquark correlations of a certain strength in- 
cluded. The reason for this is the following: The impulse ap- 
proximate photon-diquark couplings yield contributions that 
tend to fall off slower with increasing Q 2 than those of the 
quark. This is the case for both, the electric and the magnetic 
form factor. If no axialvector diquark correlations inside the 
nucleon are maintained, however, the only diquark contribu- 
tion to the electromagnetic current arises from (Jf c ) sc_sc , see 
eqs. (3.9,3.10). Although this term does provide for a substan- 
tial contribution to Ge, its respective contribution to Gm is 
of the order of 10 -3 . This reflects the fact that an on-shell 
scalar diquark would have no magnetic moment at all, and the 
small contribution to Gm may be interpreted as an off-shell ef- 
fect. Consequently, too large a ratio ji v Ge/Gm results, if only 
scalar diquarks are maintained [9]. For Set II (with weak ax- 
ialvector correlations), this effect is still visible, although the 
scalar-to-axialvector transitions already bend the ratio towards 





Set I 


Set II 


expt. 




Sc-q 


1.35 


1.33 






Ax-dq 


0.44 


0.08 




Hp 


Sc-Ax 


0.43 


0.24 






Ex-q 


0.25 


0.22 






SUM 


2.48 


1.92 


2.79 




SUM 


-1.53 


-1.35 


-1.91 


fJ>p ~\~ f^n 


isoscalar 


0.95 


0.57 


0.88 


jl>p f-i>n 


isovector 


4.01 


3.27 


4.70 



Table 4.2. Magnetic moments of proton and neutron. For 
the proton we list the following contributions separately: from 
the impulse quark-coupling with scalar nucleon amplitudes 
'Sc-q', from the axialvector diquark 'Ax-dq', from the scalar- 
axialvector transition 'Sc-Ax', and from the exchange quark 
'Ex-q'. Seagull and scalar diquark contributions are small. 



lower values at larger Q 2 . These transitions almost exclusively 
contribute to Gm, and it thus follows that the stronger axi- 
alvector correlations of Set I enhance this effect. Just as for 
the isoscalar magnetic moment, the axialvector diquark corre- 
lations of Set I tend to be somewhat too strong here again, 
however. To summarize, the ratio fi p Ge/Gm imposes an up- 
per limit on the relative importance of the axialvector corre- 
lations of estimated 30% (to the BS norm of the nucleons). 
This finding will be confirmed once more in our analysis of the 
pion-nucleon and the axial coupling constant below. 

4-2 g-nNN and qa 

Examining <7^atjv(0) which is assumed to be close to the 
physical pion-nucleon coupling at Q 2 — ~m 2 (within 10% 
by PC AC), and the axial coupling constant <?a(0), we find 
large contributions to both of these arising from the scalar- 
axialvector transitions, see table 4.3. As mentioned in the pre- 
vious section, the various diquark contributions violate the 
Goldberger-Treiman relation. Some compensations occur be- 
tween the small contributions from the axialvector diquark 
impulse-coupling and the comparatively large ones from scalar- 





Set I 

<77riV]v(0) 


9a(0) 


Set II 

<7irJViv(0) 


9a(0) 


Sc-q 


7.96 


0.76 


9.25 


0.86 


Ax-q 


0.50 


0.04 


0.10 


0.01 


Ax-dq 


1.44 


0.18 


0.34 


0.04 


Sc-Ax 


5.66 


0.39 


3.79 


0.22 


Ex-q 


1.69 


0.12 


2.70 


0.22 


SUM 


17.25 


1.49 


16.18 


1.35 


expt. 


g nNN : 13.14 ± 0.07 [25] 
13.38 ±0.12 [26] 


g a : 1.267 ±0.0035 [27] 



Table 4.3. Various contributions to <?7rjviv(0) and (?a(0), la- 
belled as in table 4.2 (with 'Ax-q' for the impulse quark- 
coupling with axialvector amplitudes) . 
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The Ratio \i p G^ 



5 Summary and Conclusions 



1.5 - 



a JLab Hall A Coll, PRL 84 (2000), 1398 [22] 
- Set I 
Set II 



eg 

O 




0.5 



0.5 1 1.5 2 2.5 3 3.5 

Q 2 [GeV 2 ] 

Fig. 4.2. The ratio (fj, p Ge)/Gm compared to the data from 
ref. [24]. 



axialvector transitions which provide for the dominant effect 
to yield g A (0) > 1. 

Summing all these contributions, the Goldberger-Treiman 
discrepancy, 



Agt 



g7riViv(0) /ir 

g A (0) M„" 



(4.1) 



amounts to 0.14 for Set I and 0.18 for Set II. The larger dis- 
crepancy for Set II (with weaker axialvector correlations) is due 
to the larger violation of the Goldberger-Treiman relation from 
the exchange quark contribution in this case. This contribution 
is dominated by the scalar amplitudes, and its Goldberger- 
Treiman violation should therefore be compensated by appro- 
priate chiral seagulls. These discrepancies, and the overesti- 
mate of the pion-nucleon coupling, indicate that axialvector 
diquarks inside nucleons are likely to represent quite subdom- 
inant correlations. 

The strong and weak radii are presented in table 4.2 and 
the corresponding form factors in figure 4.3. The axial form 
factor is experimentally known much less precisely than the 
electromagnetic form factors. In the right panel of figure 4.3 
the experimental situation is summarized by a band of dipole 
parameterizations of g A that are consistent with a wide-band 
neutrino experiment [28]. Besides the slightly too large values 
obtained for Q 2 — > which are likely to be due to the PCAC 
violations of axialvector diquarks as discussed in the previ- 
ous section (and which are thus less significant for the weaker 
axialvector correlations of Set II), our results yield quite com- 
pelling agreement with the experimental bounds. 





Set I 


Set II 


expt. 


r^NN [fm] 


0.83 


0.81 




r A [fm] 


0.82 


0.81 


0.70±0.09 



Table 4.4. Strong radius rvivjv and weak radius r A , the ex- 
perimental value of the latter is taken from [28] . 



The description of baryons as fully relativistic bound states 
of quark and glue reduces to an effective Bethe-Salpeter (BS) 
equation with quark-exchange interaction when irreducible 3- 
quark interactions are neglected and separable 2-quark (di- 
quark) correlations are assumed. By including axialvector di- 
quark correlations with non-trivial quark substructure, we 
solved the BS equations of this covariant quark-diquark model 
for nucleons and the Zi-resonance. While the A cannot be de- 
scribed without axialvector diquarks, the nucleon- A mass split- 
ting imposes an upper bound on their relative importance in- 
side nucleons, as compared to the scalar diquark correlations. 
At present, this bound seems somewhat too strong for a simul- 
taneous description of octet and decuplet baryons in a fully 
satisfactory manner. 

We furthermore extended previous studies of nucleon prop- 
erties within the covariant quark-diquark model. In this way 
we assess the influence of the axialvector diquark correlations 
with non-trivial quark substructure. Electromagnetic form fac- 
tors, the weak form factor g A and the strong form factor g W NN 
have been computed. Structures and strengths of the other- 
wise unknown axialvector diquark couplings, and of scalar- 
axialvector transitions, have thereby been obtained by resolv- 
ing the quark-loop substructure of the diquarks at vanishing 
momentum transfer (Q 2 — > 0). 

An excellent description is obtained for the electric form 
factors of both nucleons. The ratio of the proton electric and 
magnetic moment, ^jl p Ge/Gm as recently measured at TJ- 
NAF [24], is well described with axialvector diquark correla- 
tions of moderate strength. Our results clearly indicate that 
axialvector diquarks are necessary to reproduce the qualita- 
tive behavior of the experimental data for this ratio. At the 
same time, an upper bound on the relative importance of ax- 
ialvector diquarks and scalar-axialvector transitions (together 
of estimated 30% to the BS norm of the nucleons) can be in- 
ferred. 

For axialvector correlations of such a strength, the phe- 
nomenological value for the isoscalar magnetic moment of the 
nucleons is well reproduced. The isovector contribution results 
around 15% too small. While the axialvector diquarks lead 
to a considerable improvement, both magnetic moments tend 
to be around 50% too small with scalar diquarks alone [9], 
this remaining 15% mismatch in the isovector magnetic mo- 
ment seems to be due to other effects. One possibility might 
be provided by vector diquarks. While their contributions to 
the binding energy of the nucleons are expected to be negligi- 
ble, the photon couplings with vector diquarks could be strong 
enough to compensate this and thus lead to sizeable effects, in 
particular, in the magnetic moments. 

For the pion-nucleon and the axial coupling constant, we 
found moderate violations of PCAC and the Goldberger-Trei- 
man relation. For scalar diquarks alone, this is attributed to 
some violations of the (partial) conservation of the simplified 
axial current neglecting exchange quark couplings and chiral 
seagulls. Maintaining axialvector diquarks, additional PCAC 
violations can arise from missing vector diquarks which mix 
with the axialvectors under chiral transformations as pointed 
out in ref. [11]. This explains why the weaker axialvector cor- 
relations lead to better values for g A (0) and g^-iviv(O). Never- 
theless, these violations are reasonably small and occur only at 
small momentum transfers Q 2 . The axial form factor g A (Q 2 ) 
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Fig. 4.3. The pion-nucleon form factor g-n-NN(Q 2 ) and the axial form factor <?a(Q 2 ). 'Diquark diagrams' labels the sum of the 
impulse approximate axialvector contributions and scalar-axialvector transitions. The shaded region in the right panel represents 
the uncertainty in qa as determined from quasi-elastic neutrino scattering when a dipole form is fitted to both, the vector and 
the axial form factor [28]. 



is otherwise in good agreement with the experimental bounds 
from quasi-elastic neutrino scattering in [28]. 

It should be noted that qualitatively the scalar-axialvector 
transitions are of particular importance to obtain values for 
gA > 1. These transitions thus solve a problem that previ- 
ous applications of the covariant quark-diquark model shared 
with many chiral nucleon models. Their quantitative effect is 
somewhat too large, as discussed above. 

The conclusions from our present study can be summa- 
rized as follows: While selected observables, sensitive to ax- 
ialvector diquark correlations, can be improved considerably 
by their inclusion, other observables (and the nucleon-Z\ mass 
splitting) provide upper bounds on their relative importance 
as compared to scalar diquarks. These bounds confirm that 
scalar diquarks provide for the dynamically dominant 2-quark 
correlations inside nucleons. Deviations of the order of 15% re- 
main in the isovector part of the magnetic moment (too small) , 
in (?a(0) (too large) and in the Goldberger-Treiman relation. 
While these cannot be fully accounted for by including the 
axialvector diquark correlations, overall, however, the quark- 
diquark model was demonstrated to describe nucleon proper- 
ties quite successfully. 
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A Resolving Diquarks 



A.l Electromagnetic Vertices 

Here we adopt an impulse approximation to couple the photon 
directly to the quarks inside the diquarks obtaining the scalar, 
axialvector and the photon-induced scalar-axialvector diquark 
transition couplings as represented by the 3 diagrams in fig- 
ure 3.2. For on-shell diquarks these yield diquark form factors 
and, under some mild assumptions on the quark-quark inter- 
action kernel [9], current conservation followed for amplitudes 
which solve a diquark BS equation. Due to the quark-exchange 
antisymmetry of the diquark amplitudes it suffices to calculate 
one diagram for each of the three contributions, i.e., those of 
figure 3.2 in which the photon couples to the "upper" quark 
line. The color trace yields one as in the normalization inte- 
grals, eqs. (2.7) and (2.8). The traces over the diquark flavor 
matrices with the charge matrix will be included implicitly in 
those over the Dirac structures in the resolved vertices which 
read (with the minus sign for fermion loops), 

r? c = -Tr/^ x (^) S(p 2 )(-i^)S(pi) x 

x(^)s T (ps), (A.l) 

= -Trf-gfc r (^) S(P2)H 7 M )S(pi) x 

X* (^) S T ( P3 ) , (A.2) 

= -Tr/r^r X (^) S(p2)(-irnS(px) x 

X* (^) S T ( PS ) (A.3) 

= 2im g e^ ( Pd + k d ) p Q x x (A.4) 
f d 4 q g s g a P(q-Q/4)P(q + Q/A) 
J (2*)* (p? + ml)(p§ + m2)(p§ + mf) ' 
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9 r s esc /9s 


97 SC /9a 


K 


f^sa 


Set I 


0.943 


1.421 


1.01 


2.09 


Set II 


0.907 


3.342 


1.04 


2.14 





K 5 


K 5 

"7*, ax 


k 5 


K 5 


Set I 


4.53 


4.41 


3.97 


1.97 


Set II 


4.55 


4.47 


3.84 


2.13 



Table A.l. Rescaled diquark normalizations and constants Table A. 2. Strengths for pion- and pseudovector-diquark 
of photon-diquark couplings. couplings. 



The quark momenta herein are, 
Pd + k d 



Pi = 



Q , 

2 +q 



Pd + k d Q 
P2 = ^— + - + 



Pd + fed 
ps = — - A q 



(A.5) 



Even though current conservation can be maintained with these 
vertices on-shell, off-shell rjf c and Fai afi do not satisfy the 
Ward-Takahashi identities for the free propagators in 
eqs. (2.2,2.3). They can thus not directly be employed to couple 
the photon to the diquarks inside the nucleon without violat- 
ing gauge invariance. For Q = 0, however, they can be used to 
estimate the anomalous magnetic moment k of the axialvector 
diquark and the strength of the scalar-axialvector transition, 
denoted by K aa in (3.7), as follows: 

First we calculate the contributions of the scalar and axi- 
alvector diquark to the proton charge, i.e. the second diagram 
in figure 3.1 to Ge(0), with replacing the vertices r£. and r^ a/3 
given in eqs. (3.5) and (3.6) by the resolved ones, t^. and r^ a/3 
in eqs. (A.l) and (A. 2). Since the bare vertices on the other 
hand satisfy the Ward-Takahashi identities, and since current 
conservation is maintained in the calculation of the electro- 
magnetic form factors, the correct charges of both nucleons are 
guaranteed to result from the contributions to Ge(0) obtained 
with these bare vertices, rjf c and rg x a0 of eqs. (3.5) and (3.6). 
In order to reproduce the these correct contributions, we then 
adjust the values for the diquark couplings, g s and g a , to be 
used in connection with the resolved vertices of eqs. (A.l) and 
(A. 2). This yields couplings g r s esc and g a esc , slightly rescaled 
(by a factor of the order of one, see table A.l). Once these 
are fixed we can continue and calculate the contributions to 
the magnetic moment of the proton that arise from the re- 
solved axialvector and transition couplings, r^ x al3 and f 1 ^ 13 , 
respectively. These contributions determine the values of the 
constants k and K sa for the couplings in eqs. (3.6) and (3.7), 
(3.8). 

The results are given in table A.l. As can be seen, the 
values obtained for k and K 3a by this procedure are insensi- 
tive to the parameter sets for the nucleon amplitudes. In the 
calculations of observables we use k = 1.0 and K 3a = 2.1. 

A. 2 Pseudoscalar and Pseudovector Vertices 

The pion and the pseudovector current do not couple to 
the scalar diquark. Therefore, in both cases only those two 
contributions have to be computed which are obtained from 
the middle and lower diagrams in figure 3.2 with replacing the 
photon-quark vertex by the pion-quark vertex of eq. (3.14), and 
by the pseudovector-quark vertex of eq. (3.20), respectively. 
For Dirac part of the vertex describing the pion coupling to 
the axialvector diquark this yields, 

^,2 



^L = ~^ px ^(p d + kdYQ v x 



(A.6) 



/ 



A gl P(q - Q/4)P(g + Q/4) 
(27r)*(p?+mg)(pl + m?)(p§+m2) ' 



and fixes its strength (at Q 2 = 0) to n ax « 4.5, see table A. 2. 

For the effective pseudovector-axialvector diquark vertex in 
eq. (3.23) it is sufficient to consider the regular part, since its 
pion pole contribution is fully determined by eq. (A.6) already. 
The regular part reads, 



J (27T 



q g 2 gP(q-Q/4)P(q + Q/^ 
Y{pl+ml){pl + m%){pl+rnl) 



(A.7) 



-4m,e 



2 [ip\i> 



(pi +P2 +p-iY — ^757^27^17^3 



Although after the (/-integration, the terms in brackets yield 
the four independent Lorentz structures discussed in the para- 
graph above eq. (3.23), only the first term contributes to <?a(0) 
(with pi +p 2 +p 3 = (3/4) (p d + k d )+q). 

The scalar-axialvector transition induced by the pion is 
described by the vertex 



5,sa 



4*^ 
It 



d 4 q 



(P2Ps)Pi - 



,g a P(q - Q/4)P(g + Q/4) 

(P3Pl)P2 + {P1P2)PI 



(A.; 



(p\ + m%)(pl + m 2 q )(pl + m 2 q ) 



and the reverse (axialvector-scalar) transition is obtained by 
substituting Q — > —Q (or pi <-> P2) in (A. 9). The correspond- 
ing vertex for the pseudovector current reads 



r^ a = -4im q 



f d\ 

J 



g g a ga P(g - Q/4)P(g + Q/4) 
■Y{pl+rnl)(pl + w? q ){p%+ml) 

[S^ P {m q - p-Lp 2 - P2P-6 - P'APl) + 
{piP2>r +{piP3>r-{P2P3}r] ■ 



(A.9) 



The short-hand notation for a(n) (anti)symmetric product used 
herein is defined as {piP2}±" = P1P2 iPiJ 5 ^ • The reverse tran- 
sition is obtained from Q — » — Q together with an overall sign 
change in (A.9). As already mentioned in the main text, the 
term proportional to 8^ provides 99 per cent of the value for 
gA as obtained with the full vertex. It therefore clearly repre- 
sents the dominant tensor structure. 

As explained for the electromagnetic couplings of diquarks, 
we use these resolved vertices in connection with the rescaled 
couplings g r s esc and g 1 a <isc to compute g^NN and gA in the limit 
Q — > 0. In this way the otherwise unknown constants that 
occur in the (pointlike) vertices of eqs. (3.22-3.25) are deter- 
mined. 

As seen from the results in table A. 2, the values obtained 
for these effective coupling constants are only slightly depen- 
dent on the parameter set (the only exception being K^ j3a 



where the two values differ by 



lations presented in Sec. 3.2 we employ K ax = 4.5, K 5 ^ ax 



For the form factor calcu- 
4.4, 



3.9 and 



2.1. 
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